First measurement of direct /o(980) photo-production on the proton 
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We report on the results of the first measurement of exclusive /o(980) meson photoproduction 
on protons for E- t = 3.0 - 3.8 GeV and -t = 0.4 - 1.0 GeV 2 . Data were collected with the CLAS 
detector at the Thomas Jefferson National Accelerator Facility. The resonance was detected via its 
decay in the tt + tt~ channel by performing a partial wave analysis of the reaction 'yp —> p-K + ir~ . Clear 
evidence of the /o(980) meson was found in the interference between P and S waves at M n + W - ~ 1 
GeV. The S-wave differential cross section integrated in the mass range of the /o(980) was found to 
be a factor of 50 smaller than the cross section for the p meson. This is the first time the /o(980) 
meson has been measured in a photoproduction experiment. 

PACS numbers: 13.60.Le,14.40.Cs,11.80.Et 



For a long time most of our knowledge on the 
light quark meson spectrum was obtained from hadron- 
induced reactions, where typically tt, K, p or p beams 
were used, while very few studies with electromagnetic 
probes were attempted. Recently, high-intensity and 
high-quality tagged-photon beams, as the one available 
at JLab, have opened a new window into this field. On 
one hand, through vector meson dominance, the photon 
can be effectively described as a virtual vector meson. 
On the other hand, quark-hadron duality and the point- 
like-naturc of the photon coupling make it possible to 
describe photo-hadron interactions at the QCD level. 

Spectroscopy of low-lying scalar mesons is of particular 
interest. Recent advances in application of chiral effec- 
tive field theory with dispersion relations [l], S S EL HI 
led to extensive investigation of this topic. Experimen- 



tal and theoretical evidence indicates that light scalar 
mesons make a full SU(3) flavor nonet. However, the 
mass spectrum ordering of the a, K, /o(980), and ao(980) 
mesons disfavors the naive qq picture. The most natu- 
ral explanation for this multiplet with an inverted mass 
spectrum is that these mesons are diquark-antidiquark 
bound states with correct mass ordering 

@, 0, H- The 

dependence of the cross section on the momentum trans- 
fer t and resonance mass, which reflects the properties of 
the production process, might shed light on the peculiar 
structure of these mesons. For example, the authors of 
Ref. [§] suggest that a compact qq system is expected to 
be observed as a peak in the invariant mass distribution 
of the resonance decay products, while a diffuse state, 
e.g. a meson molecule, would more likely appear as a 
dip. Furthermore, the knowledge of the photoproduction 
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cross section and spin density matrix elements is relevant 
for CP and CPT violation studies via KK interferome- 
try 0. 

So far. scalar mesons have been observed in hadron- 
hadron collisions, 77 collisions and in decays of various 
mesons such as <fi, J/vP, D and B. Their cross sections are 
relatively small compared to the dominant production of 
vector mesons; however S'-wave parameters can be ex- 
tracted by performing a partial wave analysis and exploit- 
ing the interference with the dominant P- waves. The 
dominant decay mode for most of the light scalar mesons 
is the 7T7T channel. Up to now the most comprehensive 
analyses of tt + tt~ photoproduction at few GeV energies 
were performed at DESY [HI, El, SLAC 0,0] and Jef- 
ferson Lab 15[. These measurements showed the dom- 



inance of the p resonance. In the analysis of the SLAC 
data, the angular dependence was parametrized in terms 
of P-wavc alone, and no attempt was made to extract 
S-wave or higher partial waves. More recently, the HER- 
MES Collaboration investigated the interference of the 
P-wave in tt + tt~ electroproduction (with Q 2 > 3 GeV 2 ) 
with the S- and D- waves [lij . 

In this work we focus on tt + tt~ photoproduction 
at photon energies between 3.0 GeV and 3.8 GeV in 
the range of momentum transfer squared —t between 
0.4 GeV 2 and 1.0 GeV 2 and present the first analysis of 
the S'-wave photoproduction of pion pairs in the region 
of the / (980). 

The present measurement was performed using the 
CLAS detector (CEBAF Large Acceptance Spectrome- 



ter) [ljj at Jefferson Lab in experimental Hall B with 



a bremsstrahlung photon beam produced by a primary 
continuous electron beam of energy Eq = 4.0 GeV hitting 
a gold foil of 10 -4 radiation lengths. A bremsstrahlung 
tagging system with an energy resolution of 0.1% Eq was 
used to tag photons in the energy range 3.0-3.8 GeV. The 
target consisted of a 40-cm-long cylindrical cell contain- 
ing liquid hydrogen at 20.4 K. The high-intensity photon 
flux (~ 10 7 7/s) was continuously monitored during data 
taking by an e + e~ pair spectrometer located downstream 
of the target. The systematic uncertainty of the photon 
flux has been estimated to be 10%. 

Outgoing hadrons were detected and identified in 
CLAS. Charged particle trajectories were bent by a 
toroidal magnetic field (~ 0.5 T), which is generated by 
six superconducting coils. Momentum information was 
obtained via tracking through three regions of multi-wire 
drift chambers. The CLAS momentum resolution for 
charged particles is approximately 0.5-1% (a) depending 
on the kinematics. The detector geometrical acceptance 
for each positive particle in the relevant kinematic region 
is about 40%. Time-of-flight scintillators (TOF) were 
used for hadron identification. 

The interaction time of the incoming photon in the 
target was measured by detecting the outgoing particles 
in the Start Counter (ST) [I3|. It consists of a set of 24 



2.2-mm thick plastic scintillators surrounding the hydro- 
gen cell. Coincidences between the photon tagger and 
two charged particles in the CLAS detector triggered the 
recording of the events. An integrated tagged luminosity 
of ~70 pb _1 was accumulated in 50 days of running. In 
total ~20 TB of data were collected. 

The exclusive reaction jp — > p/o(980) was measured 
via the most sizable / (980) decay mode (Jo (980) — > 
tt+it- with r(7r7r)/[F(7r7r)+r(/<A')] - 75% \v^). The fi- 
nal state was selected requiring detection of both the pro- 
ton and the tt + in CLAS and reconstructing the ir~ us- 
ing the missing-mass technique. About 40M events were 
identified after all selection cuts. Calibrations of all de- 
tector components were performed, achieving a precision 
of a few MeV in the invariant di-pion mass determina- 
tion. An experimental resolution of a similar magnitude 
was evaluated from Monte Carlo simulation. 

The data analysis consisted of two main steps: 1) ex- 
traction of moments (Ylm) of the di-pion angular distri- 
butions; 2) fit of the moments with a parametrization of 
the partial waves. In the following, we briefly outline the 
procedure, referring to a more comprehensive paper (20| 
for analysis details. 

Moments Ylm(^tt) are defined as the projection of the 
production cross section on spherical harmonics with de- 
fined angular momentum L and z-component M: 



da 



dtdMmrdQn ' 

(1) 

where E~ i is the photon energy, t the invariant momen- 
tum transfer to the di-pion system squared, and M nw its 
mass. The decay angles ft^ = (O^,^) are the polar and 
azimuthal angles of the 7r + in the helicity rest frame. 

The extraction of the moments from data requires that 
the measured angular distributions are corrected by ac- 
ceptance. The CLAS acceptance and reconstruction ef- 
ficiency were evaluated with Monte Carlo simulations. 
Events were generated according to three-particle phase 
space in the same photon energy range as the experiment, 
processed by a GEANT-bascd code that included knowl- 
edge of the detector geometry and response to travers- 
ing particles, and reconstructed using the same analy- 
sis procedure that was applied to the data. Moments 
were expanded in a model-independent way in two sets 
of basis functions and, after weighting with Monte Carlo 
simulations, they were fitted to the data by maximizing 
a likelihood function. This was built on an event-by- 
event basis, to avoid binning of the experimental data. 
In the first case, the parametrization was given in terms 
of amplitudes, while in the second, moments were di- 
rectly used 2l| . In both cases the number of basis func- 
tions was limited for practical reasons. In the kinematic 
range 3.0 < P 7 < 3.8 GeV, 0.4 < -t < 1.0 GeV 2 and 
0.4 < M nn < 1.4 GeV, moments with L < 4 and \M\ < L 
were calculated as an average of results obtained with the 
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FIG. 1: Angular moments (loo) (top), (Yio) (middle) and 
(ill) (bottom) in the photon energy bin 3.4 < E-y < 3.6 GeV 
and momentum transfer 0.5 < — t < 0.6 GeV 2 . Error bars 
include the systematic uncertainty related to the photon flux 
normalization and the moment extraction procedure. The 
gray band shows the result of the fit of the moments in terms 
of partial wave amplitudes. 



two paramctrizations. 

Detailed systematic studies were performed using both 
Monte Carlo simulations and real data to ensure the va- 
lidity of the approximations and to study possible effects 
related to the basis truncation and the detector accep- 
tance. The comparison of results obtained by the differ- 
ent methods was used to estimate the systematic uncer- 
tainty related to the analysis procedure. We found that 
the variation in the moments obtained from the differ- 
ent procedures is larger than the statistical uncertainty 
and larger than other sources of systematic uncertainty, 
such as event selection cuts, detector resolution and in- 
efficiency. The final uncertainty was then obtained by 
summing in quadrature the fit uncertainty given by MI- 
NUIT, the uncertainty associated with the photon flux 
determination, and the above-mentioned uncertainty on 
the moment extraction procedure. 

The plots in Fig. [T] show the moments (loo): (^io) and 
(lii) in a selected E 1 and t bin. From Eq.[T]it is straight- 
forward to show that moment (loo) corresponds to the 
differential cross section da/dtdM„„. As expected this is 



dominated by the contribution of the p meson in the P- 
wave shown by the prominent peak at ^0.77 GeV. 
In moments (Yio) and (In), the contribution of the S- 
wave is maximum and enters via interference with the 
dominant P-wave. 

The second step of the analysis consisted of extracting 
the partial wave amplitudes from the angular moments. 
These can be expressed as bi-linear in terms of the ampli- 
tudes ai m — azrn(A, A', A 7 , E 1 , t, M n7r ) with angular mo- 
mentum I and z-projection m (in the chosen reference 
system m coincides with the helicity of the di-pion sys- 
tem): 
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where A and A' are the initial and final nucleon helic- 
ity, respectively, A 7 is the helicity of the photon, and C 
are Clebsch-Gordan coefficients. Each amplitude was ex- 
pressed as a linear combination of 7T7t amplitudes of fixed 
isospin, ai m j with / = 0, 1, 2. The number of waves was 
reduced restricting the analysis to |m| < 1, since m = 2 
waves are expected to be small in the mass range under 
investigation [2(|. The photon helicity was restricted to 
A 7 = +1 since the other amplitudes are related by par- 
ity conservation. As a consequence only three values of 
m have to be considered: m = +1, which corresponds 
to a non-helicity flip (s-channel helicity conserving) am- 
plitude, expected to be dominant, and m = 0,-1 that 
correspond to one and two units of helicity flip, respec- 
tively. In the case of the S-w&vc (I = m = 0) , only one 
amplitude is considered. The dependence on the nucleon 
helicity was simplified as follows. For a given l,ni,E 7 ,t 
set, there are four independent partial wave amplitudes 
corresponding to the four combinations of initial and final 
nucleon helicity. In general it is expected that dominant 
amplitudes require no helicity flip 



14j. On the other 



hand, we found that at least two amplitudes were neces- 
sary to reproduce the data: therefore for each l,m, with 
1 77i | < 1 we used two sets of amplitudes corresponding to 
helicity non-flip and hclicity-flip of one unit. 

For each helicity state of the target A, recoil nucleon 
A', and 7T7r system m, in a given E y and t bin, the cor- 
responding helicity amplitude ai m (s = MjL), was ex- 
pressed using a dispersion relation 22, 23|, 24( as follows: 



aimj(s) = 



Simj(s)]aim,i(s) (3) 

,Nim,l{s')p{s')aimj{s') 



da' 



where PV represents the principal value of the integral 
and p corresponds to the phase space term. In this ex- 
pression, Nim,i and Di m .i can be written in terms of 
the scattering matrix of 7T7t scattering, chosen to repro- 
duce the known phase shifts, inelasticities 0, 0], and 
the isoscalar (I = S,D), isovector (/ = P,F) and isoten- 
sor (I = S, D) amplitudes in the range 0.4 GeV < ^/s < 




the large structure at the p mass in (Yn) is due to the 
interference of the S wave with the dominant, hclicity- 



FIG. 2: Partial wave cross sections in the same kinematic bin 
as Fig. [T] The top and bottom panels show the P- and the 
5-wave, respectively. The width of the bands represents the 
uncertainty estimated as the sum in quadrature of statistical 
and systematic uncertainties as explained in the text. 



1.4 GeV. / and Si m j are matrices in the multi-channel 
space (irir, KK) relevant for the mass range considered 
in this analysis. The amplitude ai m j represents our ig- 
norance about the production process. Since discontinu- 
ities are taken into account by functions Ni m j and Di m j, 
o-imj(s) does not have singularities for s > and can 
be expanded in a polynomial function. This was chosen 
to be of second order and its coefficients are the partial 
wave analysis parameters that were extracted by the si- 
multaneous fit of the angular moments defined in Eq. [2] 
All amplitudes but the scalar-isoscalar are saturated by 
the 7T7T state. For the scalar-isoscalar amplitude, the KK 
channel was also included. In addition, to reduce sensi- 
tivity to the large energy behavior of the (irir,KK) am- 
plitudes, the real part of the integral was subtracted and 
replaced by a polynomial in s, whose coefficients were also 
fitted. The imaginary part of the integral in Eq. [3] repre- 
sents the production of long-lived (on-shell) meson pairs 
corresponding to the non-resonant part of the scattering 
process. The real part of the same integral represents 
the direct resonant production that, in the absence of 
the on-shell part, would lead to the typical Breit-Wigner 
shape. 

Partial waves a\ m up to I = 3 (F wave) were de- 
termined fitting all moments (Ylm) with L < 4 and 
< min(L,2). Results of the fit are shown as a 
gray band in Fig. [T] on top of the experimental angu- 
lar moments (Yqo), (Yio) and (Yn) in a selected _E 7 and 
t bin. As stated above, the contribution of the S-wave 
is maximum in moments (Yio) and (Yn). In particular 



non-flip wave, P m =i (A 7 = 1 



1). In moment 



(Y 10 ) the same structure is due to the interference with 
the P m =o wave, corresponding to one unit of hclicity flip 



(A, 



1 



0). A second dip near M„ n = 1 GeV is 



clearly visible and corresponds to the direct production 
of a resonance that we interpret as the /o(980). The mass 
and width of this structure are compatible with the PDG 
values (M = 980 ± 10 MeV and T = 40 - 100 MeV 0). 

It should be noted that moments of the ir + ir~ angu- 
lar distribution can be affected by baryon resonances de- 
caying to ir + p and ir~p. These contributions represent 
a background for our analysis but, having a smooth de- 
pendence on the di-pion mass, they cannot create narrow 
structures in these observables. In addition, they are ex- 
pected to be small for low moments and limited values 
of M T7r (< 1.1 GeV) that are the focus of this analysis. 

The P and S partial wave differential cross sections 
da I dtdMnn are shown in Fig. [5] As expected, the S- 
wave photoproduction is suppressed compared to the P- 
wave, which is dominated by the p meson. This can 
be explained within Regge theory because vector me- 
son production can proceed via Pomeron exchange, while 
scalars require exchange of reggeons that become sup- 
pressed as energy increases. As a test of the whole pro- 
cedure, the differential cross section da/dt for the reac- 
tion 7p — > pp — > ]37r + 7r~ was derived integrating the M W7r 
mass from 0.4 GeV to 1.2 GeV. Comparison with previ- 
ous CLAS measurements [lH and ABBHHM Collabora- 
tion data [ll[ shows good agreement, giving confidence 
in the partial wave analysis. More details can be found 
in Ref. $2$. 

The S-wave shows a clear variation in the vicinity of 
the /o(980). However, the resonance component seems to 
be embedded in a coherent background. The evidence of 
the /o(980) signal in the S'-wave is a sign that photopro- 
duction may indeed be a good tool for accessing meson 
resonances other than vector meson states [2^|. The to- 
tal S'-wave differential cross section da/dt in the region of 
the /o(980) was obtained integrating the mass in the 
range 0.98±0.04 GeV. Differential cross sections da/dt in 
Ej = 3.4 ± 0.4 GeV, for P-wave (solid dots) and S-wave 
(open circles) obtained as described above are shown in 
Fig. [3l The solid line is a prediction for the S-wave of a 
model based on Regge exchanges 

muzi. This 

was nor- 
malized to DESY K + K~ photoproduction data [28| and 
was able to reproduce the S-wave measured in the same 
channel at Daresbury [29| . The agreement of the calcula- 
tion with our data suggests that the 7r + 7r~ S-wave cross 
section extracted here is consistent with the measurement 
in the K + K~ channel. It also indicates that the present 
data can be used in phcnomcnological analyses that, ex- 
ploiting the point-like nature of photon interactions, will 
provide information about the resonance structure and 
production mechanisms. A detailed comparison between 
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FIG. 3: Partial wave differential cross sections da/dt in the 
photon energy range E 1 = 3.0 — 3.8 GeV, for the P-wave 
(solid dots) and S-wave (open circles) integrated in the 
mass range 0.4-1.2 GeV and 0.98 ± 0.04 GeV, respectively. 
The error bars include statistical and systematic uncertainties 
summed in quadrature. The line is a model prediction for the 
S- wave from Refs. SEEl. 



theory and the measured cross section will be the subject 
of future investigations. 

In summary, we measured 7r + 7r~ photoproduction in 
the photon energy range E 1 = 3.0 — 3.8 GeV and mo- 
mentum transfer range 0.4 GeV 2 < -t < 1.0 GeV 2 per- 
forming a partial wave analysis. Moments of the di-pion 
angular distribution were parametrized in terms of pro- 
duction amplitudes, expressed as bi-linear in the partial 
waves, and fitted to the experimental data. The sys- 
tematic uncertainty related to the whole procedure was 
estimated performing the analysis using different proce- 
dures and approximations. As expected, the dominant 
partial wave was found to be the one associated with the 
helicity-non-flip p(770) production. As a test, the p pho- 
toproduction cross section was extracted and found to 
be consistent with previous measurements. The interfer- 
ence between P and S waves at M Wff ~ 1 GeV clearly 
indicates the presence of the /o(980) resonance. This 
is the first time the /o(980) meson has been measured 
in a photoproduction experiment. Using a parametriza- 
tion of the individual waves based on dispersion rela- 
tions, we were able to extract the total S'-wavc differen- 
tial cross section in the mass range of this scalar meson. 
In the lowest accessible range of the momentum transfer, 
—t = 0.4 — 0.5 GeV 2 , the differential cross section was 
found to be da/dt = 0.11 ± 0.06 pb/GeV 2 , which is a 
factor of 50 smaller than the cross section for the P-wave 
integrated in the p mass range. 
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